Unconventional superconducting states on doped Shastry-Sutherland lattice 
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By using a renormalized mean- field theory, we investigate the phase diagram of t-t'-J-J' model 
on two dimensional Shastry-Sutherland lattice which are topologically equivalent to synthesized 
material SrCu2(B03)2. We find that the symmetry of superconductivity ground state depends on 
the frustration amplitude rj — t' /t and doping concentration. For weak to intermediate frustration, 
d^2„y2-wave pairing symmetry is robust in a large parameter region. Around the symmetric point 
= 1, d-wave, s-s-wave pairing as well as staggered flux may serve as ground state by varying 
the doping level. There is a first-order transition between these distinct ground states. For larger 
frustration |r;| > 1, the ground state has an s-s-wave symmetry for both hole and electron doping. 

PACS numbers: 71.10.Hf, 71.20.Li, 74.20.Mn 



I. INTRODUCTION 



Geometrically frustrated lattices have crucial im- 
pacts on the emergence of exotic electronic states in 
strongly correlated systemsi^, examples are triangular 
layered cobaltates NaxCo022,, anisotropic triangular lat- 
tice Cs2CuCl4i^ and three dimensional pyrochlore ma- 
terial K0s20g^ In particular, resonating valance bond 
(RVB) spin liquid or valence bond crystal may exist 
in frustrated quantum magnets. There is a hope that 
unconventional superconducting state may emerge upon 
doping of the frustrated magnets, it has been pointed 
out in the recent theoretical studies^'^'^'^ '^''i^^ . Recent 
discovered two dimensional synthesized frustrated ma- 
terials^ SrCu2(B03)2 is an important compound. It is 
topologically equivalent to the Shastry-Sutherlandi^iii^ 
lattice, spin-i Cu^^ lies in two-dimensional CuBOs lay- 
ers decoupled from each other by plane of Sr^+ ions, the 
antiferromagnetic exchange couplings between Cu^^ ions 
is identical to Heisenberg-hamiltonian of SS lattice and 
motivate us to investigate its doping properties. This 
lattice has been studied many years ago as a two dimen- 
sional exactly solvableAi spin model, a schematic Shastry- 
Sutherland lattice is illustrated in Fig. [TJ Let J and J' be 
the exchange couplings along the square lattice and diag- 
onal links, respectively. The production of valence-bond 
singlets on disjointed diagonal links is the exact ground 
state for J' / J > 1.477i^iii, Experiments showed that 
J'/ J — 1.574 is an optimal valuel^ for the insulator 
SrCu2(B03)2. 

There are many previous investigations on the doping 
effect of the Shastry-Sutherland lattice and various tech- 
niques have been use d^^i^°'^^i^^ . By using slave-boson 
mean-field theory^^-, the competing orders of staggered 
flux state and d-wave superconducting state are inves- 
tigated at a specific parameter regime. Similar results 
have been obtained in a recent variational Monte Carlo 
study." Based upon the analysis of t-J-V model via the 
bond-operator formulation, a number of superconducting 
states including s-wave, (s -I- id)-wave, plaquette d-wave 
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FIG. 1: Schematic structure of Shastry-Sutherland lattice. It 
includes four sublattices 1..4. The hopping integral and spin- 
spin coupling are t and J on the n.n. links (solid line) and 
t' and J' on diagonal links (dashed lines). We use a, & to 
distinguish two diagonal links with different orientations. 

are found as the ground states of the doped Shastry- 
Sutherland lattice:^ On the other hand, exact diagonal- 
ization approaches^ have been employed to study the 
ground state of finite system and no superconducting or- 
der is found to be favored on doping. 

In this paper, we apply the plain vanilla version of 
RVB theory22ii^ to study the emergence of unconven- 
tional superconductivity. We define rj = t' /t as the frus- 
tration amplitude, where t' and t are hopping integrals 
on diagonal links and square lattice links, respectively, 
and use t-t'-J-J' model to study the doping effect on the 
Shastry-Sutherland lattice. The competition among var- 
ious superconducting states will be examined for both 
hole-doping and electron-doping cases. The phase dia- 
gram is depicted as functions of rj and doping concentra- 
tion S. In particular, four distinct ground states show 
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up. We classify these states in terms of relative phase of 
mean- field pairing amplitudes. 

There are four possible ground state candidates. In 
certain limiting cases such as I77I <C 1, it is well known 
that pairing symmetry belongs to dj,2_j^2-wave with the 
superconducting order parameters on square lattice links 
Aj; = —Ay and the pairing parameters on diagonal links 
Aq = Ah = 0. Another candidate is the s-s-wave pairing 
symmetry with A^; = Aj, and A^ = Ab, while the rela- 
tive phase shift between these two distinct links is equal 
to IT. The third candidate state is staggered flux state, 
it can only be stable in negative 77 and small doping. In 
such state, the complex particle-hole mean-field param- 
eter is modulated alternatively by a staggered magnetic 
±0. The last candidate is normal metal with vanishing of 
mean-field parameters. Our calculation shows that from 
weak to intermediate frustration, d-wave state maintains 
in a large region of electron and hole dopings. Around 
the symmetric point \ri\ = 1, the symmetry of ground 
state is sensitive to the doping level since the energies of 
three distinct states, d-wave, s-s-wave pairing and stag- 
gered flux, are almost identical. For larger frustration 
\ri\ > 1, the ground state has an s-s-wave symmetry for 
both hole and electron doping. 

The rest of the paper is organized as follows. In Sec. 
II, we propose the formalism of rcnormalized mean-fleld 
theory to study the t-t'-J-J' model Hamiltonian on the 
Shastry-Sutherland lattice. In Sec. Ill, we present our 
numerical results of rcnormalized mean-field theory as 
functions of frustration and doping level, and mean-field 
phase diagram as well. Finally a summary is given in 
Sec. IV. 

II. FORMALISM 

A primitive unit cell of the Shastry-Sutherland lattice 
includes four inequivalent sites, we consider a t-t'-J-J' 
model on such lattice. The Hamiltonian can be written 
as 

H = ~ t^iH^^U,. + h.c.)P -f J^,,S, ■ S, 

(ij) 

- Ai^n,, (1) 

i 

where cj^ is to create a hole with spin a at site i, Si is 
a spin operator, fi is the chemical potential, (ij) denotes 
a square lattice or diagonal link on the lattice, tij and 
Jij stand for the hopping integrals and antiferromagnetic 
exchange couplings, respectively, tij = t and Jy = J 
on the square lattice links, while tij = t' and Jij — J' 
on the diagonal links, as shown in Fig. 1. We use t as 
an energy unit and set t/J = 3 . We choose J' / J — 
[f /tY to be consistent with the superexchange relation 
of J = it'^/U in the large Hubbard U limit. Projection 
operator.22i2i P = Y[{1 — niffiij^) removes all the doubly 

i 

occupied states. 



We define particle-particle condensate mean-field as 
well as particle-hole condensate mean-field as, 

Aij = (cj|c]j^ - cl^c'j^)o 

= (cItCjT +cI|Cjx)o, (2) 

where ()o gives expectation value corresponding to states 
without constraint of no double occupancy. Although 
the number of independent parameters in Shastry- 
Sutherland lattice is twelve, our calculation shows the 
number can be reduced to eight due to certain symme- 
try. The effect of the projection operator is taken into ac- 
count by a set of rcnormalized factorsSiiSS^ which are de- 
termined by statistical countings. Within the Gutzwiller 
approximation, the energy of physical state \ip) can be 
reduced to that of state \ipo) which is free of double 
occupancy constraint, i.e., {tplHlip) — (-0o|-ff'|V'o) — 
{ipolgtHt + gsHslipo). In homogenous case the rcnormal- 
ized factors gt = 2^/(1 -I- ^) and gs = 4/(1-1-5)^, where i5 
denotes the doping density. Thus, we have the effective 
Hamiltonian, 

■^e// = Y ~9tUj {cl^Cja + h.c) + Y 9sJijSi ■ Sj 
(y>o" {ij} 

- M^n,, (3) 

i 

and the resulting mean-field Hamiltonian can be ex- 
pressed as 

HmF = Y --^dsJtji^ijclaCja + ^tjcl^Cji, + h.c] 

-gttij {c\„Cja -I- h.c) -I- const, (4) 

with const — | Jffs [l^jj 1^ + l^y P]- We diagonalize 
the mean-field Hamiltonian (|4|) in momentum space, all 
the local order parameters and the chemical potential /i 
are self-consistently obtained for each set of frustration 
parameter 77 and doping density (5, with this procedures 
the lowest energy state can be determined. 

III. NUMERICAL RESULTS OF PHASE 
DIAGRAM AND MEAN-FIELD THEORY 

In this section, we present our numerical results 
of rcnormalized mean-field theory on the Shastry- 
Sutherland lattice. The mean-field order parameters de- 
pend on both frustration parameter 77 and doping level 5. 
In our calculations, we choose several typical frustration 
amplitude to analyze pairing symmetry for different dop- 
ing levels. Larger frustration parameter jryj corresponds 
to stronger interactions on the diagonal bonds, and the 
symmetric point \rj\ = 1 has the strongest frustration. We 
will start from the phase diagram, then provide detailed 
discussion of mean-fleld order parameters as functions of 
frustration parameter rj and doping level 5. 



3 



0.20 
0.16 
0.12 
0.08 
0.04 



0.00 



'1 

1 
1 






1 


m/ 


1 
1 
1 
I 




D 




I 


1 

s-s, 






/ 

1 






» . 




1 

1, 


s-s 



-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 




(b) - 


1 ' 1 
-- |A 1 

x,y' 













. . . . £ 



















0.00 0.06 0.12 0.18 0.00 0.06 0.12 0.18 



FIG. 2: Phase diagram of t-t'-J-J' model on the Shastry- 
Sutherland lattice as functions of doping density 5 and frus- 
tration amplitude r/. The thick solid line denotes a first order 
phase boundary while the dashed line corresponds to a second 
order transition. 



As shown in Fig. [21 there exists four distinct phases 
in the phase diagram. It is obvious that the ground 



state has 



y2-wave or d-wave symmetry in the 



limit of <C 1 at finite doping. Our results show that 
(i^2_j^2-wave state is stable in a wide parameter regions 
of — "v/OJ < T] < \/0.96 and finite doping. Previous stud- 
ies have shown the robustness of d-wave pairing against 
weak frustrations on both triangular lattices and checker- 
board lattices lii^iiSiiii. It seems that such robustness is 
universal for weakly frustrated systems. At large \r]\ > 1, 
the ground state has an s-s- wave pairing symmetry with 
Ax = Ay and ~ Ai, while the relative phase between 
Ax and Aa is tt. Recently the two families of the Fe- 
based superconductors are 1111 systems ReOFeAs with 
rare earth ions Rc^^ and the 122 systems AeFe2As2 with 
alkaline earth clement Ac^^.. An s-s-wave pairing sym- 
metry was proposed as a popular candidate for the su- 
perconducting pairing symmetry of the Fe-based super- 
conductors.^^. 

In between the above two regions, there are two non- 
superconducting states in such small parameter region 
around \ri\ = 1. The region around t] = —1 corresponds 
to staggered flux state at low doping while the normal 
metal state prevails for 77 > 1.2 at finite doping {6 > 
0.10). It is interesting to find that there is an abrupt 
change of superconducting order parameters in between 
d-wave and s-s-wave state around rj = —1 and the phase 
transition is first order. Around rj — 1 region, phase 
transition from s-s-wave to d-wave state is a weakly first- 
order transition in which parameters change continuous 
at the boundary. Moreover, the phase transition between 
staggered- flux state and d-wave state is also first-order. 
Other phase boundary corresponds to second-order. 

As we pointed out already, in the limit of weak frustra- 
tion |?7| <C 1 upon doping, the model Hamiltonian may 



FIG. 3: The magnitudes of the mean-field order parameters 
A and ^ as functions of S for (a) rj — VO.S and (b) r] — 1. 



correspond to the well-known t-J model in which the 
d-wave superconducting symmetry is the ground state. 
Our calculations are performed for various frustration 
parameter as well as doping level. In a wide range of 
parameter region, d-wave state appears to be robust as 
the ground state. In particular, our calculations show 
that d-wave state have lowest energy for positive 77 less 
than \/0.96. In Fig. 3, we present the amplitudes of 
the mean-field parameters as functions of hole density S 
for 1] = VO.S and rj = 1, respectively. As shown in Fig. 
3(a) for rj = \/0.8, a typical d-wave state is obtained and 
the parameter shows no much doping dependence. In 
the parameter region 0.96 < 77 < 1, d-wave and s-s-wave 
superconducting state are highly competing. The mean- 
field order parameters of ground state are discontinuous 
as functions of hole density 6. For better illustration, we 
take the symmetric point 77 = 1. As displayed in Fig. 
3(b), the ground state has s-s-wave symmetry at small 
doping while the d-wave state prevails for larger doping 
level. The critical doping level corresponds to 6c — 0.035. 

To reveal the competition between s-s-wave and d- 
wave states more clearly, we compare the mean-field or- 
der parameters, chemical potential as well as energy per 
site for these two states in Fig. 4 at the symmetric 
point 77 = 1. Fig. Ufa) shows parameter functions of 
d-wave, FigHJ^b) shows that for s-s-wave state in which 
I Ad I is larger than |Aa;^j,| where the subscript d denote 
the diagonal bonds. For s-s-wave, all pairing parame- 
ters change non-monotonically to zero, and then metal- 
lic state emerges smoothly. Wc plot the parameters of 
s-s-wave from 5 — 0.005, at half filling there is no self- 
consistent s-s-wave solution. Fig. ^c) shows the crossing 
of chemical potentials for those two competing states at 
the transition point 5 ~ 0.035. In Fig. ref4(d), such a 
crossing of energy per site for those two states exhibits 
itself as well. It is rather clear that a zero-temperature 
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FIG. 4: Panels (a) and (b) describe the mean-field order pa- 
rameters as functions of S for d-wave and s-s-wave state for 
r; = 1, respectively. Panels (c) and (d) correspond to the evo- 
lutions of chemical potential and energy per site as a function 
of S for two competing states. 
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FIG. 6: Amplitudes of the mean-field order parameters as 
functions of S for (a) r? = Vl.005 and (b) 77 = VI.O8. The 
inset of (a) zooms in the low doping region. 
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FIG. 5; Amplitudes of the mean-field order parameters as 
functions of 5 for rj = — \/0.95. Panels (a) and (c) correspond 
to d-wave and staggered fiux state, respectively. Panel (b) 
describes the accumulated phase of ^ for staggered flux state. 



first-order quantum phase transition may occur at the 
transition point 6 ~ 0.035. 

From Fig. [21 one can see that near 77 = —1, a 
staggered-flux state may appear in a small parameter re- 
gion. For instance, we plot the mean-field parameters 
of d-wave state and staggered-fiux state for 77 = — -\/0.95 
in Fig. [n{a),(c) respectively. Calculation shows that for 
staggered flux state is a complex value. The phase 
of is TT, and the accumulating phase of ^ij on a pane 
is independent of 77 and reduces linearly with increasing 
doping, as shown in Fig. ^h). At half- filling, it is hard 
to obtain self-consistent solution. For 7/ less than 0.04, 



the staggered flux state is stable, while in high doping 
level d-wave state has lower energy. For 77 < —1 s-s-wave 
state emerges with introduce of mobile charge, and have 
favorable energy than the staggered flux state. 

For large frustrated amplitude, the interactions on di- 
agonal bonds may play a dominate role in determination 
of superconducting pairing symmetry. When 77 takes a 
value slightly larger than 1, as the superconducting pair- 
ing symmetry may change from s-s-wave to d-wave and 
mean-field parameters varies rather smoothly. This tran- 
sition is weakly first order. In Fig. ^a), it shows that 
varies nonmonotonically to zero and ground state evolves 
from s-s-wave state to d-wave state with increasing dop- 
ing for 77 = \/l-005. Precise calculation of pairing param- 
eters shows that around the critical point A^^y ^ 0. This 
is illustrated in the inset picture of Fig. EJa), and indi- 
cates that this is a weakly first order phase transition. 
Fig. El^b) presents the mean-field parameter as functions 
of J for 77 = \/l-08. We find that a larger 77 corresponds to 
a smaller amplitude of A^^y of d-wave state. For t] > 1.4, 
amplitude of the d-wave state is vanished and metal state 
follows the s-s-wave. 

For large frustrated amplitude, s-s-wave state is the 
ground state for both positive and negative 77. Fig. [Tfb) 
takes 7] — \/2 as an illustration. By increasing doping to 
a considerable high level, both |Aa;_j,| and jA^I approach 
zero. However they do not reach zero simultaneously 
and Ax.y decreases more rapidly. It implies that super- 
conducting order parameter may exist only on diagonal 
bonds in some cases. As positive 77 becomes larger, no 
metal state will appear since the pairing parameters of 
s-s-wave state may have finite amplitude at high doping 
level. For negative larger frustration amplitude cases, 
s-s-wave state is ground state at all doping level. Am- 
plitude of Arf is larger comparing to that of correspond- 
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FIG. 7; Amplitudes of the mean-field order parameters as 
functions of 5 for (a) rj — — \/2 and (b) 77 — \f2. 



ing positive case, since the negative t' frustrated hopping 
thus enhance pairing amplitude. It indicates that su- 
perconductivity favors electron doping. This has been 
shown in Fig. [Tlja). It should be pointed out that for 
larger frustrated case, our mean field theory can not ob- 
tain the exact dimer ground state at half-filling. 



IV. SUMMARY 

We have employed the renormalized mean-field the- 
ory to study the t-t'-J-J' model on the geometrically 



frustrated Shastry-Sutherland lattice for both hole and 
electron doping cases. Our calculation shows that the 
ground state of the doped system depends on the frus- 
tration amplitude rj and doping level 5. For weak frus- 
tration ?7 << 1, d-wave state is stable in a large param- 
eter region in agreement with the case of t-J model on 
square lattice. For strong frustration 77 > 1, s-s- wave 
state dominates in a wide range of parameter region. 
This feature has also been found in the doped triangu- 
lar and checkerboard antiferromagnets. When approach- 
ing the most frustrated point 77 = 1, d-wave state com- 
petes with s-s-wave state, the phase transitions are first- 
order, the parameters change suddenly at critical point. 
Near rj — —1, staggered flux state dominates. For frus- 
trated amplitude is not very large, as doping increasing 
ground state changes from s-s-wave state to c?-wave state 
via the weakly first-order transition. Moreover, we have 
found the enhancement of superconducting order param- 
eter for negative rj because the negative t' may introduce 
frustration in kinetic energy and result in the enhance- 
ment of the pairing amplitude. Our theoretical predica- 
tions might be examined in future experiments on doped 
SrCu2(B03)2. 
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